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Directional dark matter by polar angle direct detection and application of columnar 

recombination 
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Center for Underground Physics, Institute for Basic Science (IBS), Daejeon 305-811, Korea 

We report a systematic study on the directional sensitivity of a direct dark matter detector that 
detects the polar angle of a recoiling nucleus. A weakly interacting massive particle (WIMP)-mass 
independent method is used to obtain the sensitivity of a general detector in an isothermal galactic 
dark matter halo. By using two-dimensional distributions of energy and polar angle, a detector 
without head-tail information with 6.3 times the statistics is found to achieve the same performance 
level as a full three-dimensional tracking dark matter detector. Optimum operation orientations are 
obtained for various experimental configurations, with detectors that are space- or Earth-fixed, have 
head-tail capability or not, and use energy information or not. Earth-fixed detectors are found to 
have best sensitivity when the polar axis is oriented at a 45 degree angle from the Earth’s pole. With 
background contamination that mimics the WIMP signal’s energy distribution, the performance is 
found to decrease at a rate less than the decrease of signal purity. The WIMP-mass dependence 
of the performance of a detector with various energy thresholds that uses gaseous xenon as target 
material is reported. We find that with a 5 x 10“'*®cm^ spin-independent WIMP-nucleon cross- 
section and a 30 GeV WIMP, a 770 kg-year’s exposure with a polar detector of 10 keV threshold 
can make a three sigma discovery of directional WIMPs in the isothermal galactic dark matter halo. 

For a columnar recombination detector, experimental considerations are discussed. 

PACS numbers: 95.35.+d, 29.40.-n 


I. INTRODUCTION 

Many experiments have attempted to directly detect 
Weakly Interacting Massive Particle (WIMP) dark mat¬ 
ter candidates via their elastic scattering on target nu¬ 
clei [T]. The LUX experiment has recently limited 
the spin-independent cross sections to be under 7.6 x 
10“’^® cm^ for 33 GeV WIMP mass [5]. A method that 
has commonly been used is to measure signals associated 
with the deposited energy of the nuclear recoil. However, 
background processes, such as those induced by neutrons, 
can mimic WIMP signals |3]. To overcome this problem, 
a “smoking gun” WIMP signal would be its unique di¬ 
rectional event-rate dependence [3]. When viewed from 
the Earth, the average WIMP velocity in galactic coordi¬ 
nates is from the Cygnus direction, while the directions 
of background sources are fixed in an Earth-based co¬ 
ordinate system. This effect can be seen by the direc¬ 
tions of the recoiling nuclei that are strongly correlated 
with the directions of the incoming WIMPs [5] , and pro¬ 
duce a diurnal variation of rates due the Earth’s rota¬ 
tion EE]. Prospects for a working directional detector 
have focused on low pressure gas time projection cham¬ 
bers, such as DRIFT-II [8], DMTPC [9], NEWAGE [TO], 
MIMAG [TT], and the prototype [Tl], and emulsion 
techniques m- The sensitivity of those detectors have 
been studied assuming they are either capable of recon¬ 
structing full three-dimensional (3D) tracks [5l[T4] or only 
tracks projections on a particular detector-fixed plane 
(2D detector) (SlfTS]. It is found that of order tens of 
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events are required to reject isotropy of recoil angle distri¬ 
butions, and the number of required events are one order 
of magnitude larger for a detector that cannot measure 
the sense, defined as the absolute sign {x or —x), of the 
recoil vector. However, it is still a daunting task to fully 
reconstruct sub-100 keV nuclear recoil tracks. 

In this paper, we consider a detector that can detect 
the polar angle of the recoil track with respect to a fixed 
z axis, while the azimuthal angle is not measured. In 
the following, we refer to this as a polar detector. In 
general, since reconstruction of recoiling nucleus tracks 
is not required for a polar detector, the experimental re¬ 
alization can be more reliable and feasible. One issue for 
a directional dark matter detector is whether or not it 
possesses sense recognition capability, where we call the 
former a head-tail detector and the latter an axial de¬ 
tector. Considerable technological effort is required to 
provide head-tail detection capability. Throughout the 
paper, the standard isothermal galactic halo model is 
used to model the WIMP velocity distribution. In this 
case, a head-tail polar detector with z axis aligned with 
the WIMP wind direction has the same performance as 
a full 3D detector, because in this orientation the az¬ 
imuthal distribution is completely flat and provides no 
information at all. 

One example of an axial polar detector is a stilbene 
crystal, which is an organic single crystal whose scintilla¬ 
tion efficiency depends on the nuclear recoil direction rel¬ 
ative to crystallographic axes [16] . Recently, a new tech¬ 
nique of columnar recombination (GR) has been shown to 
be capable of realizing a working axial polar detector m 
without sense detection capability. The directional sen¬ 
sitivity in a CR detector comes from the dependence of 
electron-ion recombination level on the angle between re- 
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coil track and electric field. Here the experimental signals 
are the scintillation light for the recombinations, and drift 
electrons from the surviving ionizations. 

In most previous sensitivity studies such as those re¬ 
ported in Refs. [SJ HU flMSC)] . the recoil energy is inte¬ 
grated out and only distributions with respect to angular 
variables are studied. When all the energies are inte¬ 
grated out, the shape of angular distribution does not 
depend on the WIMP mass [2T] and, thus, the distribu¬ 
tion and the analysis are greatly simplified. However, we 
note that in a real experiment, all relevant information 
are used in order to achieve the highest possible sensitiv¬ 
ity. Thus, both energy and directional information are 
used simultaneously in this work. Previous studies that 
did use both the energy and angular information [3U US] , 
treated the WIMP mass as a known parameter that was 
kept at a fixed value in the fit. While this procedure is 
simpler to implement for producing plots of cross-section 
upper limits versus WIMP mass, it is problematic for a 
general directional sensitivity calculation. In this study, 
we address this issue by treating the WIMP mass as a 
nuisance parameter in a statistical framework that uses 
the well-established profile likelihood method. 

We study signals from spin-independent WIMP- 
nucleus interactions, and assume zero background to pro¬ 
vide benchmark results. For nonzero background exper¬ 
iments, an observation of directionality still unambigu¬ 
ously leads to an observation of WIMPs, and the amount 
of directional sensitivity loss due to background is also 
studied. 

In sections that follow, we will start by calculating the 
distribution of observables, and then find the best ori¬ 
entations for various detector conhgurations using the 
standard prohle likelihood method. Later, the effects of 
background contamination are calculated for detectors 
in optimal orientations. Finally, we study the perfor¬ 
mance of detectors using xenon target with various en¬ 
ergy thresholds and WIMP masses, and estimate the re¬ 
quired detector exposure to see a three sigma directional 
signal. 


II. DISTRIBUTIONS OF OBSERVABLES 


The general double-differential derivative of the recoil 
rate per unit target mass, with respect to nuclear recoil 
energy E and solid angle H in the nuclear recoil direction 
q, can be expressed as [2^ : 
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Here /i = mM/{m -\- M) is the reduced mass of the 
WIMP-nucleus system, where m and M are the WIMP 
and target nucleus mass, respectively. E is the form fac¬ 
tor, (To is spin-independent WIMP-nucleus cross-section, 
and n = p'^/m is the number density of the WIMP, where 
we use = 0.3 GeV/cm'^. The minimum WIMP speed 


required for a recoil of energy E is u„iin- The Radon 
transform, /(umin,?), of the WIMP velocity distribution 
/(v) represents the sum of the probability densities where 
the velocity projection on direction q is equal to Vmin'- 

f(Vmin,q) = J S{v ■ q-V^in)fMd^V. ( 2 ) 


In the standard isothermal galactic halo model, the dis¬ 
tribution of WIMP velocities v relative to the target is 
Maxwellian, with an average value of ve: 


/(v) 
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In this case the Radon transform is: 
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The maximum recoil energy for a WIMP with velocity 
V is Ejnaxiv) = 2p?v'^/M, and by reversing the formula, 

we obtain Umin = \J^j^ = definition 

Eo = E^ax{vo)- We also calculate the standard total rate 
Ro when ve = 0 to be Rq = nao/M ■ J vf{v)\yj^=od^v = 
2/y^ • naoVo/M. For a a single target-nucleus type, the 
differential rate exhibits more physical meaning when it 
is expressed using Eq and Rq: 


dEdn AeEo ^ ’ 
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It is convenient to define two variables: 



In a reference frame where the reference z axis is parallel 
to the WIMP wind direction, or the Cygnus direction, 
the solid angle dependence reduces to polar angle 9 de¬ 
pendence, and the differential rate can be expressed as: 
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For simplicity, the peculiar velocity of the Sun and the 
Earth’s orbital velocity about the Sun are ignored and the 
velocity of Earth ue is taken to be equal to the velocity 
of the Local Standard of Rest 0lsr- In this case the 
magnitude of ue and vq is same, which is usually given 
the value of 220 km/s. Hence in this paper we use cce = I 
to represent the standard halo model. 

It is feasible to build a detector that is oriented at a 
fixed direction in space, where the WIMP wind direction 
and the z axis of detector form a fixed angle 9q. In the 
frame of a detector system, the polar and azimuth angles 
of recoil 9^ and (f>L provide the directional information 
and the cos 9 term in Eq. 0 can be expressed in terms 
of 9e and (j)L as cos 9 = cos 9q cos 0 e + sin 9o sin cos (j)L ■ 
Here the WIMP wind direction is assumed to lie in the 









3 



Here t is zero when the WIMP velocity lies in the plane 
made by the detector z axis and Earth’s polar axis. Since 
t distributes uniformly within 0 and 1 when all energies 
and angles are considered, the appropriate distribution 
for the Earth-fixed case, now in terms of (E,cos6*^,t), 
dEdcose dt ’ determined by replacing Oq in Eq. (j^ with 
Eq. (10) and integrating out (j)L. 

Note that here, instead of using angular information 
only as in [3 m m [Tu, we use two and three dimen¬ 
sional distributions that include energy and sidereal time 
as additional variables. 


III. STATISTICAL TEST FOR DIRECTIONAL 
SIGNAL 


FIG. 1. Schematic plot of an Earth-fixed detector shown as a 
rectangular box. 


X — z plane, without any loss in generality. The resulting 
differential rate in lab frame with respect to E, cosOl 
and is then: 


d^R 


dEd cos 9Ld(j)L 
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Here, a solid angle transform ^ „ —>• xi • 

’ ® 27z dEd COS 9 dEdcos9Ld4>L 

is applied from Eq. ([^. 

For a detector where the only directional measurable 
is the polar angle 9^, the azimuthal angle (j)L should be 
integrated out: 


The most common method that is used in particle 
physics to determine significance of an observation is the 
profile likelihood ratio test statistic method. It uses a 
hypothesis test against a null or trivial hypothesis Hq, 
where the data is assumed to correspond to the distribu¬ 
tion of an alternative or interesting hypothesis Eli. Usu¬ 
ally, a set of parameters of interest e of Hi are fixed to 
a specific value eq to obtain Hq. A likelihood ratio is 
calculated as the maximum likelihood of a null hypothe¬ 
sis divided by the maximum likelihood of an alternative 
hypothesis, with nuisance parameters 9 floating: 


A(i^o) 


L(t/ = 1^0,9) 

L{i>,9) 


( 11 ) 


The test statistic q is defined as —21nA(r'o). The p 
value is the probability to have a discrepancy larger than 
the observed one as 

pOO 

P= f{q\Ho)dqo. (12) 

Jgabs 


d^R _ d^R 

dEd cos 9L Jo dEd cos 9Ld(j)L 

For an axial detector where the sense of the recoil 
direction cannot be distinguished, the distributions 

should be folded, as ^Ed\£X\d<^b = dEd£Xd^r. + 
_ ’Em _ 'Er = I 

dEd{—cos 9 L)d{(l)L-\-'^)'’ dEd\cos9E\ dEd cos 9 l ' 

d'^R 

dEd{— cos 9l) ' 

Now we consider a common case where a detector is 
fixed in the Earth’s coordinate frame. Since the Earth 
itself is rotating, the angle 9o modulates with a period 
of one sidereal day. The WIMP wind from Cygnus is 
directed to the Earth at a constant angle e = 42° to the 
Earth’s polar axis, as shown in Fig. If the z axis of 
the detector is oriented relative to Earth’s polar axis at 
a fixed angle 0 _d, then the dependence of the angle 9o as 
a function of time t in units of sidereal days is: 

cos 00 = — cose cos 0_D — sine sin 0D cos(27rt). (10) 


The smaller the p value, the more credible that hypoth¬ 
esis Hq is not correct. For example, a p value of 0.00135 
corresponds to a Scr signal significance. In principle, the 
distribution f{q\Ho) needs to be obtained using Monte 
Carlo methods with high statistics. However, according 
to Wilk’s theorem, we can assume q follows a distribu¬ 
tion [23 for data sets in our study. Then the p value can 
be directly calculated as the probability above q°^^ for a 
distribution P with degrees of freedom equal to the 
number of fixed parameters in the numerator of Eq. (111: 
p = P{x)dx. The significance Z is related with p via 


Z = $“^(1 — p), where 4) is the cumulative function of 
the standard Gaussian. 

The likelihood function L used in this paper is the un¬ 
binned product of normalized probability density func¬ 
tion (PDF) of observed quantities for all events. The 
PDF is normalized over two or three dimensional vari¬ 
ables, such as {E,cos9l) or with the inclusion of t. A 
nonzero Xe indicates a finite average WIMP speed from 
the Cygnus direction. Thus, for studies in this paper. 
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xe as V is the sole parameter of interest, with the null 
hypothesis being xb = 0. The WIMP mass, since it is 
unknown, is a nuisance parameter and floated in the fits 
for both the numerator and denominator of Eq. ( [TT] ) . Un¬ 
like most of the studies such as |53] , where a fixed WIMP 
mass need to be inserted in order to get the sensitivity or 
limit, our method directly focuses on the directionality 
itself. 

To test the sensitivity of a measurement, the median 
(not mean) value of the test statistic qmed of a large 
ensemble of hypothesis tests is used as a measurement 
quantity of sensitivity for the work reported here. The 
ensemble size N is 1000 for all of calculations discussed 
below. The error of g^ed is 2 VNf{q^^i) ’ /(9) 

the normalized probability density function that is ob¬ 
tained from the raw distribution of q for the ensemble by 
a kernel density estimation method [26]. Here, for each 
hypothesis test, a toy Monte Carlo data set with a cer¬ 
tain number of events is generated and fitted to obtain 
the q value for this data set. 

With only one parameter of interest, the significance 
in units of a is calculated as the square root of the test 
statistic: Z = ^/q. So a Z = Scr significance corresponds 
to q = 9. In all the studies of this work, gmed is found 
to be proportional to the total event number, as will be 
shown in Fig.j^ So we use a linear proportional function 
to fit qnied, and find the number of events corresponding 
to gmed = 9 as the required size for a 3 (t discovery. Be¬ 
cause of the linearity, the sensitivity can also be scaled 
proportionally to the size of statistics when needed. 



FIG. 2. Median q values on solid markers as a function of co¬ 
sine of the angle between the detector 2 : axis and the WIMP 
wind direction for general space-fixed detectors. Head-tail po¬ 
lar and angular only detectors correspond to up and down tri¬ 
angles. Axial polar and angular only detectors correspond to 
circles and diamonds. The event number for head-tail config¬ 
uration is 200, different from 50 events for axial configuration, 
in order to fit in a single plot. 


IV. GENERAL POLAR DETECTOR 


To study the general performance of a polar detector, 
a general detector with zero energy threshold and unit 
form factor is considered. In the absence of the form fac¬ 
tor term in Eq. ([^, the energy dependence of the rate 
is expressed through the ratio E/Eq, i.e., in units of Eq. 
Since the WIMP mass dependence of the rate function 
only enters through Eq, a change of WIMP mass just 
changes Eq, the overall scale of the energy E. Thus the 
sensitivity for the general polar detector has no depen¬ 
dence on the WIMP mass. In Secs. IV m general results 
applicable to all WIMP masses are shown. 

Figure shows the median value of q for space-fixed 
detectors, with head-tail and axial configurations, for dif¬ 
ferent detector orientation angles. We can see that the 
optimum orientation for a space-fixed case is always at 
cos do = 1) where the detector’s z axis points toward the 
Cygnus. This is expected since in this case the detec¬ 
tor can gain maximum information from the measured 
cosOl- Nontrivial curves for axial detectors are discov¬ 
ered with a minimum at around cos 9q = 0.56, where the 
distribution of | cos^lI is close to flat due to the folding 
of positive and negative cos9l values. 

In the Earth-fixed case, the quantity q^ed is averaged 
over many directions, and depends on the orientation 9d. 


Figure]^ shows the dependence of qmed on the angle 9e. 
A general head-tail detector is shown to have the most 
sensitivity when the z axis is oriented parallel to the 
Earth’s polar axis (du = 0°). Interestingly, the optimal 
9e value for a general axial detector is near 9d = 45°. A 
fit using an empirical function q^ed = Q cos{a{9E — d^)) 
gives the optimum value d^ = 44.48 ± 0.33, with the 
fitted curve shown in Fig. as a dashed red line. The 
reason for d^ being slightly larger than 42 degree, is that 
it is more advantageous to lie inside than outside for the 
WIMP wind vector, with respect to the cone made by 
the detector’s z axis following the Earth’s self-rotation. 
Note that this value d^ will change if the angle of the 
WIMP wind to the Earth’s pole is a value different than 
42 degree. 

Results for the required number of events for 3cr discov¬ 
ery for four detector types, i.e., cos dp = I for space-fixed 
detectors, and 9d = 45(0)° for Earth-fixed axial (head- 
tail) detectors, are displayed in row one of Table |Tj Of 
order 10 to 20 event are required for head-tail detectors. 
These results show that a change from a space- to Earth- 
fixed basis worsens the sensitivity by factors of 3.1 and 1.9 
for axial and head-tail configurations. For an Earth-fixed 
detector, we conclude that the loss of head-tail capability 
causes the sensitivity to be a factor of 10 worse, and the 
axial one is less efficient than a full 3D one (column 2) by 
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FIG. 3. Median q values on open markers as a function of 
angle 9o between detector 2 axis and Earth’s pole, for Earth- 
fixed detectors. General head-tail polar and angular only de¬ 
tectors with 50 events correspond to up and down triangles. 
General axial polar and angular only detectors with 200 events 
correspond to circles and diamonds. For a standard axial 
xenon detector of 3 keV threshold with 200 observed events, 
four graphs for WIMP masses of 20, 30, 50 and 100 GeV are 
shown as blue squares with decreasing sizes. These graphs are 
fitted using an empirical function gmed = Qcos(a(0_D — Ol)) 
that are superimposed as red solid curves for a standard ax¬ 
ial xenon detector, and as red dashed and dotted curves for 
general axial polar and angular only detectors. 


about a factor of 19. These numbers can be compared to 
the ID readout in Ref. [53], which gives values of 8 and 
^ 10 for the two factors, by measuring the signal yield 
with a constrained background. 

As discussed in the Introduction, the optimal space- 
fixed head-tail detector behaves the same as a full 3D 
detector with its z axis aligned along the WIMP wind. 
By dividing the two numbers for the space-fixed case, we 
expect that an axial detector would have the same sensi¬ 
tivity as a full 3D detector, provided it accumulated 6.3 
times as many events. This is a much better conclusion 
than that in Ref. [53], which did not study a space-hxed 



FIG. 4. Median q values versus numbers of observed events for 
all detector types in their optimal orientations. Solid squares 
represent a standard space-fixed axial xenon detector. Other 
symbol markers are explained in Figs. and [^ The fitted 
linearly proportional lines are superimposed. Blue dashed 
lines represent axial xenon detectors configured at 30 GeV 
WIMP and 3 keV threshold. 

detector. 


V. ANGULAR ONLY DETECTOR 

Most previous estimates of the sensitivity of a direc¬ 
tional dark matter detector have only considered angular 
information only. To study detectors with polar angle in¬ 
formation only, the rate as a function of cos 0/, is obtained 
from Eq. (|^ by integrating out the energy part: 

dR r dE ndi 

d cos 9l Jo dEdcosdL' 

The same statistical procedure is performed as for the 
general polar detector case, for space- or Earth-fixed, ax¬ 
ial or head-tail configurations. The ^med dependence on 
detector orientation is shown in Figs. and In all 
four detector configurations, the shape of the cos Oq and 
Ou dependence is similar to that for the general polar 
detector case. For a space-fixed axial detector, the de¬ 
pendence on cos 00 is also not monotonic. In the Earth- 
fixed case, the red dotted curve in Fig. has a maxi¬ 
mum at 9d = 43.68 ± 0.28°. Results for the optimal 
value, cos 00 = 1 for space-fixed, 0 _d = 45° and 0° for 
Earth-fixed axial and head-tail detectors, are reported 
in row 2 of Table |T] For head-tail detectors, of order of 
10 (20) events are required to see directionality in space 
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TABLE I. The required number of point interactions for 3a discovery for all detector types. Numbers are obtained from the 
abscissa value corresponding to ^med = 9 level in Fig. for the fitted straight lines. For xenon detector types, the numbers 
are calculated by applying the ratio A^pint/Aobs from the required number of observed events. The indicated errors are purely 
statistical from finite ensemble size of 1000. 


Detector type 

space 

-fixed 

Earth-fixed 

Axial 

Head-tail 

Axial 

Head-tail 

General 

Angular Only 

Xenon, 30 GeV WIMP, 3 keV threshold 
Xenon, 30 GeV WIMP, 10 keV threshold 
Xenon, 50 GeV WIMP, 30 keV threshold 

62.30 ±0.56 

231.4 ± 2.7 
168.1 ± 1.6 

218.5 ±2.0 
1318 ± 12 

9.90 ±0.09 
12.70 ±0.13 

193.2 ±2.1 
767.0 ±9.2 

535.6 ±5.8 

702.7 ±7.6 
4087 ± 44 

19.22 ±0.19 
23.91 ±0.25 


(Earth)-fixed basis, similar to results from previous stud¬ 
ies [Himiini for 3D and 2D detectors. The sensitivity 
worsens by factors of 3.3 and 1.9 for axial and head-tail 
configurations, respectively, when changing a space-fixed 
to an Earth-fixed basis, similar to the general polar de¬ 
tector case. 

By comparing the event numbers in columns 1 and 
3 of Table i.e., for the axial configurations of general 
and angular only detectors, we can see that 3.7 and 4.0 
larger statistical samples are needed when the energy in¬ 
formation is not used. However, for head-tail detectors 
(columns 2 and 4), only 1.3 and 1.2 times the statis¬ 
tics are required. This is because in the absence of the 
powerful head-tail information, which directly relates to 
directionality, we need to rely on the correlation between 
recoil energy and direction to improve the directional sen¬ 
sitivity. Thus, the importance of combining energy and 
angular observables is established, at least for axial de¬ 
tectors. We note that the performance may be improved 
over that reported in Refs. [HI [19], if the observed energy 
is simultaneously used. 

VI. BACKGROUND EFFECT 

In nondirectional dark matter searches, the most dan¬ 
gerous backgrounds are those with the same energy- 
dependent shape as the signal, since in those cases one 
cannot distinguish the background by the energy distri¬ 
bution. In this study, we use the directional observables 
to distinguish the signal, while the directional informa¬ 
tion is diluted by the background contribution. 

Consider a background that is distributed the same 
way as the WIMP signal, while not requiring the rela¬ 
tive speed between the earth and the dark matter halo. 
By setting xe = 0 and a unit form factor in Eq. Q, 
the background rate i?BG is isotropically distributed 
with a simple exponential decaying energy dependence: 

d^Raa _ J^„-E/Eo 
dEd COS Ol 2Eq 

In a real experiment, usually we do not fully under¬ 
stand all the sources of background. Here we show that 
this is not a problem in our method. Using the same 
PDF [Eq. ([^] to fit data, regardless of the how the data 
is generated or originated as would be in a real experi¬ 
ment, a nonzero xe still unambiguously indicates a di- 



FIG. 5. Sensitivity ratio over that for pure signal case, as a 
function of signal purity defined as the ratio of signal compo¬ 
nent over all events in the data set. For each detector type, 
data points are normalized to a constant number of signal 
events. Space- and Earth-fixed polar detectors under axial 
and head-tail configurations are shown when placed in each 
one’s optimal orientation. 


rectional signal. We use the same statistical procedure, 
and only change the generation of the signal-only data 
set to a mixture of signal and background. The parame¬ 
ter of interest is still xe- Our method avoids the need for 
an external constraint to be imposed on the background 
rate, as done, for example, in Ref. [23] to obtain a better 
control of the signal component. 

The sensitivity dependence on the signal purity is 
shown in Fig. where the signal component yield is a 
constant for each detector type. All polar detectors are 
aligned in their optimal orientations. The x-axis is the 
signal purity, defined as the ratio of the input signal num¬ 
ber over total event number in the generation. In Fig. 
the sensitivity level can be seen to be almost linear with 
the signal purity. As the signal purity goes lower, the 
relative sensitivity falls by a lesser amount, especially for 
the space-hxed axial case. Thus we conclude the polar 
detector has a guaranteed performance level in the pres¬ 
ence of background. In particular, the required number 
of events for purity 0.4 only increases by about a factor 
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of 2 compared to the pure signal case. This can be com¬ 
pared to the ID case of Fig. 2 in Ref. [23], in which case 
the degradation of sensitivity is significantly worse as the 
signal purity decreases. 

Note that, other types of backgrounds, such those with 
a flat energy dependence, are easier to handle because the 
difference in the energy distribution gives extra informa¬ 
tion to distinguish them from the real signal. 


VII. STANDARD XENON DETECTOR 

Now we consider a detector with a specific material and 
a minimum energy threshold. Since the proposed colum¬ 
nar recombination detector uses gaseous xenon as target, 
a xenon Helm form factor [27] J-{E) is computed and in¬ 
serted into Eq. Q. For the minimum energy threshold, 
the LUX experiment |5] provides a good reference, where 
the trigger efficiency reaches 50% for 4 keV nuclear re¬ 
coil events. In gaseous xenon, the commonly used reflec¬ 
tive material polytetrafluoroethylene (PTFE) has a low 
refractive coefficient of about 60% for the 175 nm scin¬ 
tillation wavelength of pure xenon, compared to more 
than 90% |28j in liquid xenon. An experimental proto¬ 
type |2n] using pure xenon has put a high threshold of 
30 keV, mainly due to a low light collection efficiency of 
3%. 

A small amount of trimethylamine (TMA) mixture in 
the xenon gas, can convert the primary excitation to ion¬ 
ization of TMA molecules through the Penning transfer 
process [30], and increases the CR effect since more ion¬ 
izations will participate in the recombination process m- 
At the same time, the scintillation light is wavelength- 
shifted to the TMA emission band centered at near-UV 
(300 nm), where the reflectance of PTFE increases to 
90% [31]. Thus, it should be possible to achieve an en¬ 
ergy threshold much lower than 30 keV. 

In this study, we consider three minimum energy 
threshold values of 3, 10 and 30 keV. Because of the 
rapidly falling shape of the form factor, the event rate 
will be suppressed by a factor more than 10“® for recoil 
energies greater than 100 keV. Thus, the maximum of 
the detection energy range is chosen to be 100 keV. 

This standard xenon detector can be viewed as a 
columnar recombination detector with perfect resolution. 
So we only study the axial configuration. First, in the 
Earth-fixed case, the sensitivity dependence on 9d is 
studied for four typical WIMP mass values of 20, 30, 
50 and 100 GeV. Figure]^ shows the results for a 3 keV 
threshold, with fitted curves superimposed. The fitted 
optimal values for Od are 45.08±0.22,45.21±0.30,44.12± 
0.49,44.30 ± 0.79 degree for 20, 30, 50 and 100 GeV 
WIMP mass. These are all nearly equal; a weighted av¬ 
erage value is 44.97 ± 0.17. In the following, the optimal 
directions 9d = 45° and cos0o = 1 are used for Earth- 
and space-fixed standard xenon detectors. 

Because of the detection threshold, the required num¬ 
ber of observed events for detectors with different thresh- 
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FIG. 6. Median q values as a function of WIMP mass for axial 
xenon detectors, normalized to 200 events of point interaction. 
From the top to bottom, the dashed, dotted and dot-dashed 
lines corresponds to the 3, 10 and 30 keV energy threshold. 
The best orientations on space-fixed (cos So = 1, main plot) 
and Earth-fixed {9d = 45°, inset plot) basis are used. 


olds cannot be directly translated to detector perfor¬ 
mance, since in this case the total rate depends on the 
threshold. The number of observed events, Aobs, is the 
integral of the differential rate within the detection range 
including the form factor for a given total target mass 
and measurement time T: 


Nohs = RMtT = 


r (PR 

'det.range dEdEl 


MtT. (14) 


To make fair comparisons of the sensitivities, we intro¬ 
duce a quantity called the number of point interactions, 
A^pint, which is the number of WIMP-nucleus interactions 
with no threshold and no form factor effects: 

poo p 

= l.SOiRqMxT = 2.608/pTT ■ no’QVoMj’T/M. 

(15) 


The amount of required Apint directly relates to detec¬ 
tor performance. For general detectors discussed in pre¬ 
vious sections, Aobs = Apint. For a single target, the 
spin-independent WIMP-nucleus cross-section ctq is re¬ 
lated to the WIMP-nucleon cross-section as (Tq = 
where ^iw-n is the WIMP-nucleon re¬ 
duced mass. For xenon, Apint = 1-70 x 10"^^anMxT/m ■ 
GeV • (cm^ • kg • year)“^. 

There will be some WIMP-mass dependence of the sen¬ 
sitivity since the form factor and the threshold depends 
on absolute energy values that do not scale with Eq. The 
WIMP mass is still a floating nuisance parameter in the 
fitting, while being fixed at various input values in event 
generation. Figurej^shows the raw sensitivity as a func¬ 
tion of WIMP mass, for space- and Earth-fixed detectors 
at their optimal orientations. Each point in the plot is 










normalized to A^pmt = 200, so that the sensitivities can 
be compared for the same detector exposure. 

In general, the higher the threshold, the lower the sen¬ 
sitivity. The decrease of the sensitivity for low WIMP 
masses is due to the effect of the threshold that elim¬ 
inates a higher fraction of events the lower the WIMP 
mass. The sensitivity decrease for high WIMP masses is 
due to the form factor, which suppresses the high energy 
recoils. The most sensitive WIMP-mass value increases 
as the energy threshold increases, and is around 30 GeV 
for a 10 keV threshold and 50 GeV for a 30 keV threshold. 
The inset of Fig.j^shows that the performance of Earth- 
fixed detectors follows the same behavior as space-fixed 
case with respect to WIMP-mass and energy-threshold 
dependence, except for an overall sensitivity that is lower 
by a factor of about three, which is consistent with the 


general study in Sec. IV 


We can also compare the decrease in sensitivity to the 
decrease in the observed event rate due to the energy 
threshold cut. For example, for a 30 GeV WIMP, an en¬ 
ergy threshold change from 3 keV to 10 keV, decreases 
the sensitivity by only 23%, compared to the 71% de¬ 
crease in the total number of observed events. This is 
due to the fact that most directional information comes 
from the high recoil energy events, so that a cut on low 
energy events has minor impact. For a detector with 
high energy threshold of 30 keV, a WIMP mass in the 
intermediate range around 50 GeV is most promising, 
with the sensitivity that is 3.2 times below that of a 3 
keV-threshold detector for the same range. 

Last three rows of Table |T] show the required point 
interaction numbers for a 3(7 discovery under several 
WIMP-mass and energy-threshold combinations. In a 
typical case of a 30 GeV WIMP mass and a 10 keV en¬ 
ergy threshold, the numbers are 219 and 703, for space- 
and Earth-fixed xenon polar detectors. These correspond 
to 770 and 2480 kg-year exposures for a 5 x 10““^® cm^ 
spin-independent WIMP-nucleon cross-section. 


VIII. DISCUSSION 

A columnar recombination detector using high pres¬ 
sure xenon-TMA gas mixture can measure the angle be¬ 
tween the nuclear recoil track and electric drift field [HI 
[32], without sense recognition capability. Thus it is an 
axial polar detector. The columnar effect benefits from 
a high Penning efficiency as the percentage of the xenon 
excitations that fall back to free electrons by ionizing the 
TMA molecules. For electron and gamma energy deposit, 
the Penning efficiency is estimated to be around 10% and 
20% in Ref. [33] and Ref. [34], and also indirectly mea¬ 
sured to be 10%-15% in high pressure xenon gas with 
TMA mixture [35j . In addition, the primary scintilla¬ 
tion from xenon is observed to be absorbed by the TMA 
mixture [35]. Nevertheless, the columnar recombination 
effect for the ionization channel has been seen for the al¬ 
pha particle events as shown in Fig. 5 of Ref. [33], where, 


as a first approximation, the collected charge is a linear 
function of cos^ 9^. Further studies are needed to obtain 
a high Penning efficiency for the high initial ionization 
density produced by nuclear recoils. 

To date, no experimental results on the properties of 
nuclear recoils in a columnar recombination material is 
available. Once a reliable detector response is available, a 
detector resolution can be straightforwardly included in 
our model to obtain the sensitivity. Despite this limita¬ 
tion, we have tested a naive resolution smearing method 
for nuclear recoils based on the following considerations: 
(a) the recombination fraction of the electron-ion pairs 
ranges from 0.8 for recoils that are parallel(0L = 0°) to 
0.4 for those that are perpendicular {9l = 90°) to the 
field direction; (b) an ionization work function of 76 eV 
for the electron-ion pairs; and (c) 10% and 50% efficiency 
for the recombination (scintillation) and ionization chan¬ 
nels with Gaussian resolutions calculated from a Fano 
factor m of 0.14 [35] for xenon gas. The resulting sensi¬ 
tivity is only worsened by 7% relative to the zero resolu¬ 
tion case. This can be understood as that the direction¬ 
ality is derived from an overall anisotropic phenomenon 
and not from a localized spot. Reference [20] also noted 
that variations in angular resolution do not make much 
difference in the dark matter detection sensitivity. Thus 
we do not expect a significant change to our sensitivity 
results by detector resolution effects. 


IX. CONCLUSION 

The performance of a directional dark matter detector 
with polar angle detection is studied for various config¬ 
urations. A WIMP-mass independent method is used to 
obtain the sensitivity of a general detector. In addition, 
a detector with xenon as target material is studied. We 
infer that: 

• Both axial and head-tail polar detectors have the 
highest sensitivity when the z axis is aligned with 
WIMP wind. However, the dependence of sensi¬ 
tivity to detector orientation is not monotonic for 
an axial detector. To obtain optimal performance 
when rotating with the Earth, the z axis should be 
oriented at 45 degree to Earth’s pole for an axial 
detector, while it should be aligned with Earth’s 
polar axis for a head-tail detector. 

• A head-tail polar detector can detect directionality 
with of order 10 or 20 events on a space- or Earth- 
fixed basis. In the absence of sense detection ca¬ 
pability, an order of magnitude more statistics is 
needed. 

• Without using energy information simultaneously, 
the required statistics would be a factor of 3.7 (4.0) 
times higher for a space (Earth)-fixed axial detec¬ 
tor. This conclusion will be useful for detector 
types in which partial directional information is 
available, such as a 2D planar detector. 
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• A general axial polar detector with 6.3 times the 
statistics has the same performance as a general 
full 3D tracking detector. However, in experimen¬ 
tal practice, the target mass for a full 3D detec¬ 
tor is limited because of diffusion effects, and ac¬ 
complishing millimeter tracking is extremely chal¬ 
lenging. On the other hand, a detector with polar 
angle sensitivity without head-tail discrimination 
requirement can use straightforward experimental 
techniques. In addition, it can be made in large 
volume with high-density gas, with a target mass 
that can be orders of magnitude larger than a con¬ 
ventional full 3D detector. We conclude that it is 
of great advantage to explore the directional dark 
matter detection technique using a senseless polar 
angle detection apparatus. 


• A space-fixed detector is generally found to be 3 
and 2 times more sensitive than an Earth-fixed de¬ 
tector, for axial and head-tail configurations. This 
ratio is an important factor when comparing the 
additional cost for a space fixed detector, since it 
has to rotate all the time with respect to the Earth. 
In a space-fixed detector, to distinguish between 
the WIMP signal from the galaxy coordinate and 
possible anisotropic background originating from 


the detector frame, manually reversing the detector 
z axis direction for half of the measurement period 
should be useful. 

• In the presence of contamination by an isotropic 
background that mimics the energy shape of the 
signal, the decrease of directional sensitivity occurs 
at a rate that is less than the decrease of signal pu¬ 
rity. Thus a polar detector is robust against back¬ 
ground. The space-fixed axial configuration is least 
affected by such backgrounds. 

• For an axial xenon polar detector, the decrease of 
sensitivity due to the effect of the energy threshold 
is much smaller than the decrease of the total event 
rate. A 770 or 2480 kg-year exposure can reach a 3a 
directional signal on a space- or Earth-fixed basis 
with 10 keV energy threshold, for a 5 x 10“^® cm^ 
WIMP-nucleon cross-section and a 30 GeV WIMP 
mass. 
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